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Abstract 

Formation of chaos in the parametric dependent system of interacting oscillators for the both 
classical and quantum cases has been investigated. Domain in which classical motion is chaotic is 
defined. It has been shown that for certain values of the parameters from this domain, form of the 
classical power spectrum is in a good agreement with the quantum band profile. Local density of 
states is calculated. The range in which application of perturbation theory is correct has been defined. 

Pacs number 05.45. -a 

1 Introduction. Problem statement. 

At last time the investigation of parametric dependent Hamiltonians like H(Q, P, A) is a subject of great 
interest. Here (Q, P) is a set of canonical coordinates and impulses, and A is a parameter describing 
connection of the system with the external field. Interest to such systems is motivated by mesoscopic 
physics [1]. Most of the works devoted to the parametric dependent systems are dealing with the in- 
vestigation of following situation: For A = Hamiltonian is exactly integrable. For A > Hamiltonian 
H(Q, P, A) becomes non-integrable and for some defined values of A = Ao, solutions of classical equations 
generated by the Hamiltonian H (Q, P, Ao) display chaos. After this, the little variation of the parameter 
SX = A — Ao is of interest. Using the basis |<j>„) in which H(Q,P,Xq) is diagonal, for small values of <5A, 
Hamiltonian H(Q, P, A) can be written in the form of sum of two matrix 

H = E + SXB. 

Here Eq is the diagonal matrix, and B is the banded matrix, off-diagonal elements of which are random 
numbers [2], By way of numerical diagonalization it is possible to define eigenvectors |$„(Ao + <5A)), 
|$„(Ao)) and eigenvalues E n (\ + 5X),E(X ) corresponding to the Hamiltonians H (Q, P, X),H (Q, P, Ao) 
appropriately. The main interest at this is connected with the dependence of parametric kernel 

P{n\m) = |($„(Ao + ^A)|$ m (A ))| 2 



from the parameter <5A. The value of P(r) — P(n\n + r) averaged over the n can be considered as the 
local density of states. Descriptions of foregoing methods of solving of problems of quantum chaos are 
given in [3-5] . Goal of this paper is investigation of chaos in the system of interacting oscillators when 
both of the Hamiltonians Hq(Q, P, Xq), H{Q, P, A) display chaos. Situation when the Hamiltonian Hq is 
integrable has been studied by us earlier (see papers [6,7]). 
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2 Classical consideration. Model Hamiltonian. 



Investigation of the following Hamiltonian is of interest for us 
H = H Q + XV, 

p2 p2 p2 -i 
tj r I , r 2 , r 3 , ^ f 2 , 2 , 2\ 

H o = Y + y + Y + 2™ +?2+93)> 
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Here and below we work with the dimensionless quantities. Canonical equations for the Hamiltonian 
are of form 



(2) 



'7i 


= Pi, 




Pi 


= -qi 


- ZMiqha 
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= P2, 
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= -92 


- 2\q2qfql 


<k 


= P3, 




Pa 


= -Q3 


- 2Xq 3 qfql 



Solutions for the set of equations in case of A ^ can be obtained only numerically. As is shown by 
numerical analyze, solutions of the set of equations J2J f° r the values of the parameters A = 3.0, H = 4.2 
are essentially chaotic (see Fig.l). For further calculations we have to obtain classical power spectrum 
with the purpose of its further comparison with the band profile. Let us construct correlation function 
for V(qi q_2 qa) 



dH 

dx 



G(r) = ((V(t + r)-V) (V(t)-V)) 



(3) 



where (...) means averaging over the time, V is the averaged values of V. We have to make Fourier 
transformation of the expression (J3J to get classical power spectrum 



C(w) 



e -n^G{r)dr. 



In the practical calculation of Q we use the method equivalent of (@J 



C(f) = 



T 



T/2 



T/2 



'V(t)dt 



(4) 



(5) 



here x £ [0, T] is time interval in which solution of the set of equations J2J) is obtained, / = To get 
the power spectrum one needs to use discrete form of expression (JSJ 



C{k) 



TA 



21 



n— — l 



(6) 



i2Tvk 

)e "4 



k = fNA, N is a number of dots dividing the interval T, A = I is the 



here X(fc) = Ejto 1 ^0% j - *° " xx^^x wx ^v^ i& »<*x ^ , — - ^ 

averaging parameter [8,9]. Form of the classical power spectrum for the different values of the parameter 
I is given on the Figures Fig. 2, Fig. 3. These data are confirming the fact that solutions of J2J) are chaotic. 
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3 Quantum mechanical consideration 



Our subsequent goal is the comparison of the band profile | B nm \ 2 considered as the function of u> = " E "^" E "' 
with the classical band profile. Let us set about to the quantum-mechanical analyze of the Hamiltonian 
(Q. Eigenvectors of the basis, in which Hamiltonian H(P, Q, 0) is diagonal can be found easily and 
are of the form [10] 

^ („,„„„„\ - ^i^E^M^^ e-hiil+il+il) ( 7 ) 

^nin 2 n z \<J1<1213) ~ , , , j j j — T e j { ' ) 

^/2 ni + n2 + n3 niln 2 ln 3 lTr 3 

where H n (q) denotes the Hermite polynomials [11]. Eigenvalue energy is given by the formula 

E n = (n+l), (8) 

where 

n = n\ + n 2 + n 3 . (9) 
Condition JjJJl implies existence of the ( n+1 H"+ 2 ) f ^ degeneration for the given n. 

Diagonalization of the Hamiltonian (1) can be done only numerically. On the Figures 4, and 5 dependence 
of the energy eigenvalues E n (X) from the parameter A are plotted. According to Fig. 4 for the A ^ 
removal of degeneration takes place. In case of further increasing of the parameter A up to the values 
characterizing the domain of chaotic motion A = 3.0, as easy to see from Fig. 5, repulsion of the energy 
terms happens. This fact is of great importance since the repulsion of the energy terms is a sine of 
emerging of quantum chaos [12]. For the defining of the band profile at first we need to find eigenvectors 
of the Hamiltonian H(P, Q, Xq). This problem is possible to resolve by way of numerical diagonalization 
of the Hamiltonian H(P, Q, Aq), at the values of Xq = 3.0 in the basis of functions Q 

\(H + X V)\* > (io) 

Since the matrix elements are satisfying the condition 

(n \V(qiq 2 q 3 )\ m) = {n 1 n 2 n 2 \V{qiq 2 q 3 )\ mim 2 m 2 ) = (m \ql\ mi) (n 2 \q 2 \ m 2 ) (n 3 |gf| m 3 ) , 
taking into account that 

q nm = -j^Wn8 n -\,m + Vn + lS n+ i >m ), (11) 
and using the condition (q 2 ) nm = J2k QnkQkm we get 

(9 2 )«m = \(VnVn- IS n -2, m + (2n + l)5 n ,m + Vn + lVn + 25 1l+2 , m ). (12) 

In l|ll|) and (|12|1 5„ >m stands for Kronecker delta function. Substituting l|llf> and (|12|) into (|10|) after 
straightforward calculations one can obtain 

1 



nm 2 n 3 \V\ mim 2 m 3 ) = -(vWv^i = Wm-2,mi + ( 2 "i + l)S numi + Vm + Wni + 2S ni+2>mi ) x 



X (y/n2~ \Jn 2 - l(>„ 2 _2, m2 + (2^2 + l)^n 2 ,m 2 + V"2 + 1a/«2 + 25„ 2+2 ,ma) X 
X (y/n3~Vn 3 - l(5„ 3 _2,m 3 + ( 2 «3 + l)5n 3 ,m 3 + \/™3 + Wn 3 + 25 n3+2 , m3 ), 

3 

(nin 2 n 3 \H \mim 2 m 3 ) = (n 1 + n 2 + n 3 + (13) 

After forming from the matrix elements l|13|) matrix and making diagonalization one can obtain eigen- 
values spectrum E n /x ) and basis of eigenvectors $ n (Ao). Calculating matrix elements of V(qi, q 2l q 3 ) in 
the basis $„(Ao) we can get the band profile 

Bin = ^-\(*M\V(qi, 92,93)1 <f m (A ))| 2 (14) 
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Here Aq = ah d is the mean level spacing, and d is the system's dimension. Coefficient a is determined 
from the condition a = jjrm > where D(E) = |^ is the density of states, and SN is the number of states 
in the energy interval SE = E(P, Q, Ao + SX) — E(P, Q, Xo). Density of states in it's turn is possible to 
determine from the numerical diagonalization of the Hamiltonian H(P, Q,X$) , or from the integral 

D{E) = — / dPdQ. (15) 

{2irh) a dE J H (P,Q,\)=E 

Here the integral must be taken over the phase space enclosed by energy surface (see Fig. 6) In our 
calculations we are interested in the domain in which classical motion is essentially chaotic: E(P, Q, X = 
3) = 4.2; SX = 0.1; SE = 0.1; % = 0.1; a ~ 7.2. Taking smaller values of Ti leads to the necessity of 
diagonalization of matrix larger then our. In our case {% — 0.1) matrix dimension is 455 x 455. Band 
profile B^ m as a function of uj = is plotted on the Fig. 7. According to Fig. 7 band profile is in a 

good qualitative agreement with the classical power spectrum. Better quantitative agreement is possible 
to achieve for smaller values of fi [1] (since the quasi-classical approximation implies K — ► 0). Another 
quantity we are interested in is the parametric kernel 



P(r) = |($„ +r (Ao + <5A)|$ n (A ))| 2 , (16) 

where the averaging is done over the several n states. Evidently (|16f) determines correlation between the 
wave functions. The parametric kernel l|16|l frequently is named as the Local Density of States (LDS) 
[3]. We are interested in P(r), since with its help, is possible to determine the boundaries in which 
perturbation theory is valid. For that one has to compare parametric kernel (|16|) with the results for 
LDS obtained by means of perturbation theory [5] 

St 2 \B I 2 

where the parameter T(Sx) for given Sx must be defined from the normalization condition 

^P prt (r) = l. (18) 

r 

Results of numerical calculations for the values Sx = 0.1 and Sx = 0.2 are presented on the Figures Fig. 8 
and Fig. 9. On basis of obtained results we conclude that boundary in which perturbation theory for our 
system is correct is limited by the condition Sx < 1. 
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Fig.l Phase trajectory projection on the plane (pi,qi) obtained for the values of the parameters: 
E = 4.2, A = 3.0. It is easy to see that trajectory is chaotic. 
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Fig. 2 Classical Power Spectrum plotted for the values of the parameters: E = 4.2, A = 3.0, 1 = 10 




Fig. 3 Classical Power Spectrum plotted for the values of the parameters: E = 4.2, A = 3.0, I = 3 
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Fig. 6 Equipotential surface constraint by the condition H{P, Q, A) = E = 4, 2 
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Fig. 7 Quantum band profile B^ m as a function of u> = En ^ Em for the values of parameters 
E = 4.2, A = 3.0, SX = 0, 1. Evidently B\ m is in a good qualitative agreement with the Classical Power 

Spectrum (see Fig. 2 ) 
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Fig. 8 Quantum Profile P(r) (Local Density of States) and First-Order Perturbative Profile P pr t(r) for 

the values of the parameter A = 0, 1. 
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Fig. 9 Quantum Profile P(r) (Local Density of States) and First-Order Perturbative Profile Pprtix) for 
the values of A = 0, 2. Comparing Fig. 8 with Fig. 9 one could say that agreement between P(r) and is 

better if A < 0, 1 . 
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